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ABSTRACT 

We present a new technique for constraining the topology of the universe. The method 
exploits the existence of correlations in the phases of the spherical harmonic coefficients 
of the CMB temperature pattern associated with matched pairs of circles seen in the 
sky in universes with non-trivial topology. The method is computationally faster than 
all other statistics developed to hunt for these matched circles. We applied the method 
to a range of simulations with topologies of various forms and on different scales. A 
characteristic form of phase correlation is found in the simulations. We also applied 
the method to preliminary CMB maps derived from WMAP, but the separation of 
topological effects from, e.g., foregrounds is not straightforward. 

Key v^rords: cosmic microwave background - cosmology: theory - methods: statistical 



1 INTRODUCTION 

The Ist-year Wilkinson Microwave Anisotropy Probe 
(WMAP) observations seem in good accord with the emerg- 
ing standard cosmological model; a flat A-dominated uni- 
verse seeded by scale-invariant adiabatic Gaussian fluctua- 
tions (Spergel et al. 2003). However, at large scales there is 
an unexpected loss of CMB anisotropy power that was also 
seen in the COBE-DMR data (Bennett et al. 2003a). One 
possible explanation is that we inhabit a universe that has 
a non-trivial topology. That is to say, our universe is in fact 
multi-connected and has a finite volume. If so, then power 
on scales exceeding the fundamental cell size will be sup- 
pressed. Our space may not be large enough to support long- 
wavelength fluctuations. A number of authors have tried to 
restrict the topology of the universe using the WMAP data 
(Luminet et al. 2003; Cornish et al. 2004; de Oliveira-Costa 
et al. 2004; Roukema et al. 2004; Uzan et al. 2004; Weeks 
et al. 2004). In particular, Luminet et al. (2003) have shown 
that the Poincare dodecahedral space {Qo ~ 1.013) accounts 
for the measured power in the quadrupole and octopole bet- 
ter than an infinite (simply-connected) fiat universe. The 
angular power spectrum, however, is not an effective way 
to characterise the peculiar form of the anisotropy manifest 
in small universes of this type (Levin, Scannapieco & Silk 
1998). 

Whether we can determine the topology of the universe 
depends on its volume. If the universe is small enough, we 
should be able to see right around it, since photons can 
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travel across the whole universe. If so, then we may be able 
to identify ghost images of the same object in different di- 
rections in the sky or recognise the topology from signatures 
in clustering statistics. Luminet & Roukema (1999) provide 
a review of these methods. However, these methods are hin- 
dered by the need to identify good standard candles; objects 
that can be traced through different eras. This is not a prob- 
lem when looking at the CMB for signatures of the topology. 
The CMB photons originate from the same epoch and from a 
very thin shell, the last scattering surface (LSS), which is the 
same when viewed from either side. If the physical dimen- 
sions of the universe are less than the diameter of the LSS 
then the sphere self intersects; the loci of self-intersections 
are circles (Cornish, Spergel & Starkman 1998). We should 
therefore be able to match patterns of hot and cold spots 
around circles. This result holds no matter how complex the 
topology. A further advantage of using the CMB as an in- 
dicator is that the last scattering surface marks the edge of 
the visible universe, making it a powerful tool for looking 
for non-trivial topology. 

In this paper, we introduce a new method to search 
for evidence of a finite universe in all-sky CMB maps. Our 
method is based on the properties of the phases of the (com- 
plex) coefficients obtained from a spherical harmonic expan- 
sion of all-sky maps, specifically we look for phase correla- 
tions arising from matched pairs of circles. In the next sec- 
tion we briefly introduce some basic ideas in topology and 
discuss the simulations with non-trivial topologies used later 
on. In Section 3 we sketch the procedure developed to de- 
tect phase correlations in the CMB. In Section 4 we discuss 
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the results of applying the method to the simulations and 
WMAP data. We draw our conclusions in Section 5. 



2 NON-TRIVIAL TOPOLOGY 

Topology deals with connectivity. To a topologist, a cof- 
fee cup and a doughnut arc equivalent, while a coffee cup 
and a bowl are distinct. In the cosmological setting, we are 
concerned with the connectivity of space: a manifold is de- 
scribed as simply-connected if any closed path can be con- 
tracted to a point. The possibility that our universe maybe 
multi-connected was first suggested by Schwarzschild (1900), 
but has often been overlooked in favour of the simplicity 
offered by trivial topological spaces. Indeed there is a com- 
mon misconception that the shape of the universe can be 
determined from Einstein's field equations, given a value of 
Qo- General Relativity only specifies the local curvature of 
space-time, so nothing in it forbids a (global) non-trivial 
topology. The Cosmological Principle merely restricts us to 
manifolds with constant curvature. Indeed, any detection of 
non-trivial topology could determine the sign of the spatial 
curvature since its characteristics (lifter for the cases of Eu- 
clidean, spherical and hyperbolic manifolds. 

A 3-dimensional manifold can be described by identi- 
fying faces of a fundamental cell/polyhedron. The WMAP 
data suggests that the Universe is very nearly spatially fiat, 
with no = 1.02 ± 0.02 (Bennett et al. 2003a). We therefore 
restrict our attention to non-trivial topologies with a fiat ge- 
ometry. The technique should nevertheless be applicable to 
spherical and hyperbolic spaces where matched circles are 
expected. 

In this study we use 6 compact orientable fiat mod- 
els that can be constructed either by identifying the sides 
of a parallelepiped or a hexagonal prism. The simplest one 
(Model 1) is the hypertorus (3-torus), which is obtained 
from a parallelepiped with pairs of opposite faces glued to- 
gether. This manifold is built out of a parallelepiped by iden- 
tifying x—>-x + h,y—>-y + b and z z + c. The next three 
manifolds are variations of the hypertorus involving identi- 
fications on opposite faces of a twisted parallelepiped. One 
of these has opposite faces identified with one pair rotated 
through TT (the twist torus; Model 2). The next identifies 
opposite faces with one rotated by -it/2 (the 7r/2 twisted 
torus; Model 3). The last of these three is obtained by pro- 
ceeding with the following identifications: (x, y, z)—^{x + h, 
—y, —z) corresponding to translation along x and rotation 
around x by tt; next {x, y, z)—>-{—x, y + b, —(z + c)) corre- 
sponding to translation along y and z followed by rotation 
around y by tt; and finally (x, y, z)~^{~{x+h), — (y+6), z+c) 
translation along x, y and z followed by rotation around z 
by TT (the triple twist torus; Model 4). Two other topologies 
are built out of a hexagon by identifying the three pairs of 
opposite sides, while in the z direction the faces are rotated 
relatively to each other by 27r/3 {2n/3 hexagon; Model 5) 
or by 7r/3 (7r/3 hexagon; Model 6). 

The simulations presented here are based on those 
in Rocha et al. (2004): they have Oa = and a 
Harrison Zel'dovich Gaussian spectrum. The simulations re- 
flect topologies with equal-sided physical dimensions that 
lead to a suppression of the quadrupole with respect to the 
high-order modes (Weeks et al. 2004). We define a dimen- 



sionless topological scale j of a simulation as the ratio be- 
tween the width of the fundamental cell and the horizon 
size. 

These simulations include only those temperature fiuc- 
tuations generated by the Sachs- Wolfe effect. Matched cir- 
cles occur because we are looking at the same point on LSS 
from different directions. In order for this to be true, the 
temperature fluctuations need to be generated at the LSS 
which is true in this case, but the SW effect only dominates 
at large scales. We therefore limit our investigation to I < 20. 
Even at these scales, however, there are three main factors 
that could confuse the statistic : (i) velocity perturbations 
generating anisotropics at the LSS: (ii) the integrated Sachs- 
Wolfe effect due to time varying potential wells crossed by 
the photons; and (iii) Galactic foreground contamination of 
CMB data. 



3 TESTING FOR PHASE CORRELATIONS 

The temperature fluctuations in the CMB at any point in 
the celestial sphere can be expressed in spherical harmonics 
as 

oo m — + i 

A{e, <t>)=J2Yl ai,mYl,m{e, 0), (1) 
1=1 m=—l 

where the o;,m, arc complex coefficients that can be written 

ffli.m = |ai,m| CXp['t(/);,,„]. (2) 

In orthodox cosmologies the temperature fiuctuations con- 
stitute a statistically homogeneous and isotropic Gaussian 
random field. In this case the phases 4>i^m are independent 
and uniformly random on the interval [0,27r] and the vari- 
ance of aim depends only upon I. Departures from ortho- 
doxy lead to differences in behaviour of the oj.m. For ex- 
ample, in anisotropic Gaussian fiuctuations the variance of 
the ai^rn depends on m (Ferreira & Magueijo 1997; Inoue & 
Sugiyama 2003). Coles et al. (2004) developed a diagnostic 
of departures from the standard assumption that involves 
the randomness of phcises. The main component of the tech- 
nique involved using Kuiper's statistic from an available set 
of phase angles. First the phases are sorted into ascending 
order, to give the set {^i, . . . , On}. Each angle 9i is divided 
by 27r to give a set of variables Xi, where i = 1 . . . n. Prom 
the set of Xi we derive two values yS'+ and S~ where 

S'+=max(i-Xi,--X2,...,l-X„| (3) 
Inn J 

and 

S'-=max(xi,X2--,...,X„-^^|. (4) 
L n n ) 

Kuiper's statistic, V , is then defined as 

V = {St + S-) -(^1 + 0.1^^+^^^ . (5) 

The form of V is chosen so that it is approximately inde- 
pendent of sample size for large n. Anomalously large values 
of V indicate a distribution that is more clumped than a 
uniformly random distribution, while low values mean that 
angles are more regular. In order to remove any artifact from 
the choice of coordinate firame that the o;,™, are measured 
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in, V is calculated for randomly rotated coordinate systems. 
The rotated coefficients are found by employing the Wigner 
D function 

0(,m = ^ ai,m' D^m,m' ("> /^i 7) ■ (6) 
m' 

where the Euler angles a, (5, 7 define the magnitude of suc- 
cessive rotations about the coordinate axes. Consequently, 
a distribution of V is obtained from one set of phases. This 
distribution is compared, using a test, with distributions 
obtained in a similar manner from 1,000 Monte Carlo (MC) 
sets of (pi^rn drawn from a uniformly random distributiou. 
The probability 'P(x^) oi obtaining a lower value of is 
calculated from the fraction of Xmc that are less than the 
obtained from the phases. A set of angles is assumed to 
be non- random if P(x^) > 0.95. 

The random-phase hypothesis can be further scruti- 
nised by investigating subsets of the phases. Those subsets 
should also be uniformly random on the interval [0,27r]. In 
this paper, we look at two subsets: (i) the phase differences 
for fixed values of m (0(+i,m — 4'i,m) and (ii) the phase dif- 
ferences for fixed values of I — 0!,m) (even Z-modes 
only). Both subsets arc of particular interest since Cornish et 
al. (2004) indicate that matched circles arc associated with 
phase correlations. In their paper, the significance level for 
detection of matched circles in the WMAP data is calcu- 
lated from 'scrambled' versions of the data. In the scram- 
bled versions, phase correlations are removed by randomly 
exchanging the ai^m at fixed I. Also, as previously mentioned 
multi-connectedness breaks the global isotropy and some- 
times the global homogeneity of the Universe (except the 
case of the projective space). This will induce correlations 
between the a;,m of different I and m. For instance, due to 
the symmetries of the hypertorus case, (ai^rniii m') ¥= ^ =^ 
m — m' = mod (2) and I — I' = mod (2) (Riazuelo et al. 
2004). These non-zero off-diagonal terms will induce phase 
correlations as detected in our study. 

Overall, for a given temperature map, we obtain a value 
of V{x^) for each mode; 18 values for subset (i) and 10 
for subset (ii). To improve the presentation of the results, 
we combine V{x^) for each subset in two ways. First, we 
count the number of modes with Vijx^) > 0.95 and find 
the mode that displays the highest value. Secondly, we per- 
form a Kolmogorov-Smirnov (K-S) test on the distribution 
of ■P(x^) over the modes. If the phases are random, the set 
of 'P(x^) should be uniform in the interval [0, 1]. To quan- 
tify the significance of the K-S value obtained, 10,000 MC 
sets of P{x^) arc generated and the probability of obtaining 
a lower K-S value is calculated. The reason for doing both 
these things is that one would expect one in twenty modes 
to yield a value of V{x^) > 0.95. The second approach gives 
a (very conservative) idea of the significance of the whole 
set of modes rather than each individual one. 



4 RESULTS AND DISCUSSION 

The simulations were generated in HEALPix^ format 
(Gorski, Hivon & Wandelt 1998) with a resolution param- 
eter iVside=32. The a;,^ were derived using the 'cmafast' 

t http://www.eso.org/science/healpix/ 



routine in the HEALPix package. V was binned from 0- 2.75 
and 0-2.5 for subsets (i) and (ii), respectively. Subset (i) re- 
quired 10,000 rotations in order to obtain stable results. To 
analyse one realisation took 18 minutes on 1,400 MHz CPU 
desktop. On the other hand, 3,000 rotations produced stable 
results for subset (ii), resulting in each analysis taking 2| 
minutes on the same desktop. 

The six flat models listed in Section 2 were studied with 
j=0.5. In order to see the effect j had on the results, the 
hypertorus was explored in more detail. A further six simu- 
lations with j <2 were scrutinised. Also, six simulations in 
which matched circles are not anticipated {j > 2) were stud- 
ied. The results for subsets (i) and (ii) are shown in Tables 
1 and 2 respectively. The tables show the number of modes 
with 'P(x^) ^ 0.95, the most non random mode and the 
K-S fractional probabilities. For subset (i), the method was 
applied to five realisations (sets of rotations) of each simu- 
lation. The average values obtained from these realisations 
are shown in the tables. For subset (ii), only one realisation 
was necessary in order to obtain consistent results. 

If the phases are random, the number of modes with 
■P(X^) > 0.95 should be 0.9 for subset (i) and 0.5 for subset 
(ii) . From Table 1 it is clear that we are finding correlations 
when scanning across fixed m for most of the simulations (we 
shall refer to these as m-correlations) . The m-correlations 
are less obvious in terms of the K-S probabilities. However, 
this test is more general: it does not search specifically for 
large values of ^(x )) so it should be interpreted as a very 
crude measure of departure from uniformity. The average 
count for all the simulations displayed in the Table 1 is 2.5. 
The average value is significantly larger than the expected 
value of 0.9. This contrasts with the average value of 0.7, 
from Table 2, that is only slightly higher than the expected 
value. From the results, it is evident that no particular mode 
can be chosen to look for signs of non-randomness in the 
phases. This is unsurprising, as the correlations would man- 
ifest themselves across many modes, whose nature would 
depend on the location and size of the circles with respect 
to the observer. 

The spread in the number count from one realisation to 
another is quite small (roughly ±1 for larger values) and an 
exact value for each simulation can be obtained by increas- 
ing the number of realisations. To see if there is any worth 
in doing this, we looked at five further simulations for each 
model with ;;'=0.5, the results of which are shown in Table 
3. The standard deviation is very large, with a count of 1.0 
being consistent with all the models. This indicates that it 
would be very difficult to use the method to distinguish be- 
tween models or to determine the exact topological scale, 
you could merely indicate the most probable case. These 
results confirm that the hypertorus and the 7r/3 hexagon 
models, at this scale, do not display any m-correlations in 
the phases. However, the triple twist torus is detectable from 
m-correlations. This fact is less clear in Table 1, again indi- 
cating that the test results vary from simulation to simula- 
tion. 

Looking at results for the hypertorus displayed in Tar 
ble 1, it is hard to perceive any trend between the number 
count and the topological scale of the simulation. The m- 
correlations are seen both with values of j < 2 and j>2. In- 
tuitively, one would expect phase correlations when matched 
circles are present {j < 2). In such cases, the number of re- 
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peated patterns (circles in the sky) increases with decreas- 
ing j, and hence we would expect an increase in number 
count with decreasing j. This is not seen in the results, al- 
though, this may be meisked by the number count varying 
from simulation to simulation for fixed values of j. Applying 
the statistic to further simulations may reveal such a trend. 

On the other hand, the high number count seen at j > 2 
is less easy to explain. Unearthing signs of non-trivial topol- 
ogy beyond the horizon size, suggests the technique is po- 
tentially a powerful tool. Diagnostics that hunt for circles in 
the sky are limited to a maximum value of j=2. Our detec- 
tions beyond the horizon size arc by no moans serendipitous. 
Phillips & Kogut (2004) compute the covariance matrix of 
the a;,m for the hypertorus at various topological scales. 
They find the off-diagonal terms, that incorporate the phase 
information, remain prominent even when the width of the 
fundamental cell is greater than the diameter of the LSS. 

Apart from the hypertorus, each of the topologies ad- 
dressed are not only anisotropic, but also inhomogeneous. 
The question therefore arises whether the results are affected 
greatly by changes in the position of the observer. Wo gener- 
ated simulations with the observer's position shifted within 
the fundamental cell. For models 2 to 6, five simulations with 
j=0.5 were generated with the observer position shifted by 
(0.1a;,0.1a;,0.1a;) where x is the width of the fundamental 
cell. The rmmbcr count for each model was found to be con- 
sistent with those of the centrally located observer displayed 
in Table 3. 

A positive detection of m-correlations in CMB data 
would therefore be indicative of the universe having a non- 
trivial topology. In order to seek evidence of these correla- 
tions in the CMB data, we turned to four WMAP-derived 
maps. The temperature maps were all constructed in a man- 
ner that minimises foreground contamination and detector 
noise, leaving a pure CMB signal. The ultimate goal of these 
maps is to build an accurate image of the LSS that cap- 
tures the detailed morphology. Following the release of the 
WMAP 1 yr data, the WMAP team (Bennett ct al. 2003b), 
Tegmark, de Oliveira-Costa & Hamilton (2003), Nciselsky et 
al. (2003) and Eriksen et al. (2004a) have released CMB-only 
sky maps (see papers for details) . We shall refer to these as 
the ILC, TOH, Naselsky and Eriksen maps respectively. It 
is worth pointing out that even though we limit ourselves to 
full sky maps in this paper, the method can be adapted to 
cut maps following a suitable adjustment of the MC Simula^ 
tions. Nevertheless, the current WMAP data is preliminary, 
so we reserve this treatment for future data releases. Again, 
the results are shown in Tables 1 and 2. The method was 
applied to five realisations of each CMB map for subset (i) . 

In three of the four maps, evidence was found for the 
sort of phase correlations seen in the simulations. However, 
all the maps, bar the Naselsky map, display Z-correlations 
that were discussed in Coles ct al. (2004) which can be ex- 
plained, at least partially, by foreground contamination. If 
no m-correlations were seen then this would suggest there 
is no evidence for non-trivial topology. However, a positive 
detection leaves open the possibility of a non-trivial topol- 
ogy. 



Map 


Count (Mode) 


K-S Probability 


Model 1 j=0.5 


0.6 


0.68 


Model 2 j=0.5 


1.2 (17) 


0.35 


Model 3 i=0.5 


4.6 (4,5) 


0.99 


Model 4 i=0.5 


0.6 


0.22 


Model 5 i=0.5 


1.8 (14) 


0.68 


Model 6 i=0.5 


0.2 


0.42 


Model 1 j=OA 


5.4 (1,4,9,12) 


0.99 


Model 1 j=0.6 


0.4 


0.63 


Model 1 j=0.8 


0.8 


0.35 


Model 1 j=0.9 


0.0 


0.56 


Model 1 j=1.0 


0.4 


0.47 


MnHel 1 i— 1 fi 


Q a CS Q 101 


0.77 


Model 1 j=2.0 


1.4 (16) 


0.43 


Model 1 j=4.0 


0.0 


0.66 


Model 1 j=5.0 


2.8 (12) 


0.48 


Model 1 j=6.0 


0.0 


0.62 


iViOQCl 1 J — o.U 


U.o 


U.zo 


Model 1 i=10.0 


2.1 (12) 


0.88 


ILC 


1.2 (4) 


0.55 


TOH 


0.6 


0.28 


Naselsky 


2.2 (2,4) 


0.50 


Eriksen 


2.0 (1,4) 


0.33 



Table 1. Phase correlations when scanning across fixed m. The 
column labelled Count shows the average number of modes ex- 
ceeding 95 per cent significance. The column labelled Mode gives 
the mode with greatest significance. The last column shows a 
rough measure of significance for all modes obtained using a K-S 
test as described in the text. 18 modes studied in total. 



Map 


Count (Mode) 


K-S Probability 


Model 1 j=0.5 





0.17 


Model 2 i=0.5 


1 


0.66 


Model 3 i=0.5 





0.52 


Model 4 j=0.5 


1 


0.45 


Model 5 j=0.5 





0.46 


Model 6 i=0.5 


2 (2,12) 


0.13 


Model 1 j=0.4 


2 (8,20) 


0.32 


Model 1 j=0.6 





0.21 


Model 1 j=0.8 





0.68 


Model 1 j=0.9 


1 


0.97 


Model 1 i=1.0 


1 


0.08 


Model 1 i=1.6 





0.09 


Model 1 j=2.0 


1 


0.49 


Model 1 j=4.0 


1 


0.47 


Model 1 j=5.0 


1 


0.16 


Model 1 j=6.0 


1 


0.54 


Model 1 j=8.0 


1 


0.73 


Model 1 j=10.0 





0.43 


ILC 


2 (14,16) 


0.50 


TOH 


3 (6,14,16) 


0.77 


Naselsky 


1 (14) 


0.81 


Eriksen 


2 (6,16) 


0.80 



Table 2. Phase correlations when scanning across fixed I. 
Columns are as in the previous table. 10 even modes studied in 
total. 
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Map Count {Mode) 

Model 1 0.5 ± 0.2 

Model 2 1.0 ± 0.9 

Model 3 3.4 ± 3.4 

Model 4 1.0 ± 0.7 

Model .5 1.3 ± 0.8 

Model 6 0.8 ± 0.5 



Table 3. The variation in number count from simulation to sim- 
ulation. The columns show the average number of modes (along 
with the variation) and the most significant mode. We evaluated 
6 simulations with j=0.5 for each model. 

5 CONCLUSIONS 

In this paper, we have presented a new method for con- 
straining the topology of the universe. The method relies on 
utiMsing phase correlations associated with matched pairs 
of circles in the CMB sky. Wc applied the method to vari- 
ous simulations with non-trivial topologies. The method ap- 
pears to detect non-trivial topologies beyond the horizon 
size. However, the method fails to estimate the scale of fea- 
tures it detects and is not good at discriminating between 
different models. This can potentially bo overcome by study- 
ing more simulations as the results were shown to vary from 
simulation to simulation. The method is simple, computer 
tionally fast and does deliver a clear signature: a positive 
detection of these m-correlations is clear evidence for non- 
trivial topology. 

With this in mind, the method was applied to four 
CMB-only sky maps; we found evidence for m-correlations 
in three of them. However, it would be premature to con- 
clude that there is evidence for non-trivial topology in any 
of the available CMB temperature map. The WMAP data 
is preliminary and has already been shown to have a num- 
ber of unusual properties that are not yet fully understood 
(Chiang ot al. 2003; Dineen & Coles 2004; Erikscn ot al. 
2004b). Indeed, Eriksen et al. (2004a) have pointed out that 
the techniques used in producing these maps result in a poor 
reconstruction of the cosmological phases which may inter- 
fere with the possibility of detecting correlations of the type 
discussed hero. Nevertheless, with improved data, wc believe 
the method will be a useful tool in determining the shape 
of the universe. Moreover, our method is based on only the 
simplest possible measure of randomness in the phase dis- 
tribution. More sophisticated combinations may allow us to 
improve the method substantially, perhaps unearthing fur- 
ther signs of non-trivial topology beyond the horizon size. 
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